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Abstract. We present new determinations of the cosmic expansion history from 
red-envelope galaxies. We have obtained for this purpose high-quality spectra 
with the Keck-LRIS spectrograph of red-envelope galaxies in 24 galaxy clusters 
in the redshift range 0.2 < z < 1.0. We complement these Keck spectra with 
high-quality, publicly available archival spectra from the SPICES and VVDS 
surveys. We improve over our previous expansion history measurements in Simon 
et al. (2005) by providing two new determinations of the expansion history: 
H{z) = 97±62 km sec"! Mpc"! at z ~ 0.5 and H{z) = 90±40 km sec"! Mpc"! 
at z ~ 0.8. We discuss the uncertainty in the expansion history determination that 
arises from uncertainties in the synthetic stellar-population models. We then use 
these new measurements in concert with cosmic-microwave-background (CMB) 
measurements to constrain cosmological parameters, with a special emphasis 
on dark-energy parameters and constraints to the curvature. In particular, we 
demonstrate the usefulness of direct H{z) measurements by constraining the dark- 
energy equation of state parameterized by wq and Wa and allowing for arbitrary 
curvature. Further, we also constrain, using only CMB and H{z) data, the number 
of relativistic degrees of freedom to be 4 ±0.5 and their total mass to be < 0.2 eV, 
both at la. 



1. Introduction 

Direct supernova measurements of the deceleration parameter pQ, as well as indirect 
measurements based upon a combination of results from the cosmic microwave 
background (CMB) [5], large-scale structure (LSS) [3111], and the Hubble constant 
[5] indicate that the expansion is accelerating. This suggests either that gravity on 
the largest scales is described by some theory other than general relativity and/or 
that the Universe is filled with some sort of negative-pressure "dark energy" that 
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drives the accelerated expansion [Hj; either way, it requires new physics beyond 
general relativity and the standard model of particle physics. These observations have 
garnered considerable theoretical attention as well as observational and experimental 
efforts to learn more about the new physics coming into play. 

The simplest possibility is to extend Einstein's equation with a cosmological 
constant, or equivalently, to hypothesize a fluid with an equation-of-state parameter 
w = p/ p — —1 (with p and p the pressure and energy density, respectively). However, 
it may well be that the cosmological "constant" actually evolves with time, in which 
case w —1, and there are a variety of theoretical reasons [7] to believe that this 
is the case. Precise measurement of w{z) (with, in general, a parameterized redshift 
dependence) or, equivalently, the cosmic expansion history, has thus become a central 
goal of physical cosmology [H [10] . 

Among the techniques to determine the cosmic expansion history are supernova 
searches, baryon acoustic oscillations (BAO) [TTJ HH [TS] , weak lensing [T7], and 
galaxy clusters [T^. These techniques all have different strengths, and they all also 
suffer from a different set of weaknesses. As argued in the ESO/ESA and Dark Energy 
Task Force reports O [10] , robust conclusions about the cosmic expansion history will 
likely require several avenues to allow for cross checks. There may also still be room 
for other ideas for determining the expansion history. 

A common weakness of supernova searches, BAO (at least the angular clustering), 
weak lensing, and cluster-based measurements is that they rely largely on integrated 
quantities. For example supernovae probe the luminosity distance, 



The other probes rely on similar quantities, which depend on an integral of the 
expansion history, to determine the expansion history, rather than the expansion 
history itself. The purpose of the differential-age technique [18] is to circumvent 
this limitation by measuring directly the integrand, dt/dz, or in other words, the 
change in the age of the Universe as a function of redshift. This can be achieved 
by measuring ages of galaxies with respect to a fiducial model, thus circumventing 
the need to compute absolute ages. From Galactic globular clusters age-dating we 
know that relative ages are much more accurately determined than absolute ages 
(e.g., Refs. [m [2Ql [21] ) . A preliminary analysis, with archival data, has already been 
carried out [22l I23j. and the results applied to constrain dark-energy theories (e.g.. 



The challenge with the differential-age measurement is to find a population 
of standard clocks and accurately date them. There is now growing observational 
evidence that the most massive galaxies contain the oldest stellar populations up to 
redshifts ofz - 1-2 |25l[26l[2Il[28l[29]. Refs. [28] and |30] have shown that the 
most massive galaxies have less than 1% of their present stellar mass formed at z < 1. 
Ref. [29] shows that star formation in massive systems in high-density regions — i.e., 
galaxy clusters — ceased by redshift z ^ 3 and Ref. [31] shows that massive systems, 
those with stellar masses > 5 x 10^^ M©, have finished their star- formation activity 
by z ~ 2. There is thus considerable empirical evidence for a population of galaxies, 
harbored in the highest-density regions of galaxy clusters, that has formed its stellar 
population at high redshift, z > 2, and that since that thime this population has been 
evolving passively, without further episodes of star formation. These galaxies trace 
the "red envelope," and are the oldest objects in the Universe at every redshift. The 
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Refs. [24]). 
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differential ages of these galaxies should thus be a good indicator for the rate of change 
of the age of the Universe as a function of redshift. 

Here we report the first results of an observational campaign to obtain high signal- 
to-noise Keck spectra of red-envelope galaxies in a number of galaxy clusters at redshift 
z < 1.0. In particular, we report spectroscopy for galaxies in 24 rich galaxy clusters. 
These new data extend earlier results from archival field galaxy data at z < 0.5 p2l[23] 
to the redshift range 0.5 < z < 1 crucial for determining w. They also allow us to 
probe massive galaxies in clusters, which are older, and thus provide a much more 
efficient avenue to determine the age at any given redshift. 

Our paper is organized as follows. We first estimate in ij2] the statistical 
uncertainties in the age determination that arise from noisy spectra and finite sample 
sizes. The analysis demonstrates that the age and metallicity of galaxies can be 
disentangled from the measured spectra given high enough signal-to-noise, wide 
enough spectral coverage and accurate modeling, and it helps identify the wavelength 
ranges required to optimize the age determination. These arguments motivate the 
Keck observations we have made. We discuss the stellar-population models in ij3] and 
estimate the uncertainties in the age determination that arise from uncertainties in the 
stellar-population models. Section |4] discusses the observations and the galaxy spectra 
used in our analysis. These include both new spectra obtained at Keck (described in 
more detail and catalogued in Paper II, [S]), as well as galaxy spectra from archival 
data. Section [5] presents results for our measurement of the age-redshift relation, 
or equivalently, for the expansion rate H(z) as a function of z. In f}6l we combine 
our measurements of H{z) with CMB data to derive constraints to cosmological 
parameters. We conclude in §7. 

2. Statistical Uncertainties in the Age Reconstruction 

We first quantify the statistical error, arising from noisy spectra and from finite sample 
sizes, in the age determination. To do so we assume that our theoretical stellar 
models are an accurate representation of reality and then recover the galaxy age and 
metallicity from a spectrum over a given wavelength range and with a given level of 
photon noise. We have thus generated a series of mock spectra using the Chariot & 
Bruzual (2007) (CB08, Chariot private communication) stellar-population models with 
random noise for different values of signal-to-noise and for different spectral coverages. 
In Fig. [T]we explore how well we recover the age of the input spectrum for a signal-to- 
noise S/N = 10 per resolution element of 3 A. The age range is divided into 60 equally 
spaced bins in age from 1 to 14 Gyr, and 25 equally spaced bins in metallicity from from 
0.2 to 2.5 Zq. In each bin we generate 100 Monte Carlo realizations of the spectrum 
with the chosen sigma-to-noise ratio. For each realization, parametrized by a given 
age and metallicity, we maximize the likelihood in the age and metallicity parameter 
space and determine the age marginalizing over the metallicity. The metallicity is 
allowed to change from 0.2 to 2.5 Zq. The (rest frame) wavelength range used goes 
from 2500 to 7500 A. 

By comparing the long-dashed and the solid grey areas, which are the la region 
of the Monte-Carlo simulated errors, we see that the recovery is excellent and without 
bias. Even for a S/N of only 10, the age-metallicity degeneracy is lifted when blue 
light from (rest frame) 2500 A is observed. 
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Figure 1. For sufficiently high S/N data over sufficient wavelength range, the 
age-metallicity degeneracy is lifted. Input mock spectra are from the library by 
CB08 to which we have added random noise with S/N = 10 per resolution element 
of 3 A. The thin solid line is the input age. When we ffi the same CB08 spectra 
to these mock spectra with both the age and metallicity as free parameters, we 
find a very good recovery of the age, with a random dispersion smaller than a few 
percent. The age-metallciity degeneracy is removed from spectra with blue-light 
coverage (2500 - 7500 A) and modest S/N. 



3. Stellar-Population Uncertainties 

We next investigate the uncertainty that arises in the age determination from 
uncertainties in the stellar-population models. To do so, we take the mock spectra 
from fj2] with the CB08 models and then recover their age (marginalizing over the 
metallicity) with a different set of stellar-population models. Here, we use the Bruzual 
& Chariot 2003 [31] models which rely on a different set of stellar interior models 
(Geneva instead of Padova used by CB08) and different stellar atmosphere models. It 
should be noted that this exercise will generally over-estimate the uncertainty, because 
the models may have treatments of different phases of stellar evolution that are not up 
to date, or because some may be missing phases of stellar evolution [33j. As models 
get updated, they will hopefully converge toward a more unified picture [34] , 

Figs. [1H3] show the comparison between the two different models for different 
metallicities and for different wavelength coverages. As in Fig. [l] the shaded grey 
regions correspond to the la region from 100 Monte-Carlo realizations. Fig. [2] uses a 
(rest frame) spectral range from 2500 to 7500 A while Fig. [3] drops the blue light and 
covers only the range from (rest frame) 4500 to 7500 A. In both cases the S/N is 20 
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Figure 2. Age recovery for two different stellar-population models, CB08 and 
BC03, which rely on different input stellar physics (see text). The (rest frame) 
wavelength range of the fitted spectrum is 2500 to 7500 A. Note the good age 
recovery despite different input physics, especially at higher metallicities where 
the absorption lines are stronger and help lift the age-metallicity degeneracy. 

per resolution element of 3 A. Inspection of Fig. [2] shows that there is good agreement 
in the age recovery for all ages and metallicities despite the different stellar input 
physics in the models. The recovery is especially good at higher metallicities, which 
corresponds to the typical metallicity of the old, massive galaxies found in clusters. 
This is expected because the stronger metallic absorption lines provide a better lever- 
arm to determine the metallicity and thus break the age-metallicity relation. Note 
that the new CB08 models like Maraston [35] and SPEED [371 ISS] , apply a correction 
to properly model the giant branches. 

3.1. Importance of blue wavelength coverage 

It is apparent from comparison of Figs. O and [3] that the wavelength coverage of 
the observations is crucial, as the agreement degrades quickly with observations that 
do not include light bluewards of the 4000 A break. The blue end is particularly 
important, as line-blanketing at the blue end allows a more effective breaking of the 
degeneracy between the age and metallicity. It is this consideration that motivates the 
use of the Low Resolution Imaging Spectrometer (LRIS, [13]) on the Keck I telescope, 
which is senstive down to 3200 A. 

Because the method relies only on measuring relative ages, it is useful to illustrate 
the sensitivity of the observed spectrum to this parameter. In Fig. [4] we show the 




Figure 3. Same as Fig. [2] but for a (rest frame) wavelength range of 4500 to 
7500 A. Note that it becomes increasingly diffcult to recover the age when blue 
light in the spectrum is not measured. 



relative flux of spectra for different ages with respect to a 2 Gyr old solar metallicity 
population. All spectra have been arbitrarily normalised to unity at 7500 A. Even 
by eye the differences at a single wavelength are notable, above 10%. The figure also 
illustrates how the greatest sensitivity comes from observations at the blue end of the 
spectra, « 3000 A. 

4. Galaxy Spectra Sample 
4.I. Keck observations 

Motivated by the above considerations, we have started a program to obtain spectra 
of bright cluster galaxies with the Keck/LRIS instrument. We targeted rich galaxy 
clusters in order to obtain a sample as large as possible of red galaxies at redshift 
z < 1. Most of the clusters are well-known, rich X-ray clusters from a variety of 
samples such as the Abell catalog f38j, the ROSAT Cluster Survey (RCS; [39]) and 
the Massive Cluster Survey (MACS; i40j). In the redshift range 0.5 < z < 1, fewer rich 
X-ray clusters are known, so we targeted (and confirmed) two of the richest Spitzer 
mid- infrared selected cluster candidates from the IRAC Shallow Cluster Survey |41j . 
Many of the targeted clusters were also observed in the near-infrared cluster survey 
of [42], which provided a valuable and consistent astrometric resource for slitmask 
designs. The list of clusters targeted is presented in Table 13.11 and an example of 
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Figure 4. Relative flux of theoretical single stellar populations (SSPs) in the 
wavelength range 2500 to 7500 A for four different ages: 3, 5, 7, 9 Gyr, from top 
to bottom as a ratio of the 2 Gyr spectrum. The plot illustrates how differential 
ages are measured and the essential role of blue wavelength coverage for this 
experiment. 



observed spectra is shown in Fig. [5l 

The key goal of this program is to provide high-signal-to-noise, wide-wavelength- 
coverage spectroscopy of a large number of early-type galaxies at moderate redshifts. 
These spectra are then modeled to derive the ages of the galaxy stellar populations. 
Since rest-frame UV light probes the youngest and most massive stars in a galaxy, 
blue sensitivity is crucial for this experiment. Of all the optical spectrographs on 
8 — 10m class telescopes currently available, LRIS on the Keck I telescope is unique in 
being the only dual-beam spectrograph, thus providing sensitive observations across 
the entire optical window (A ^ 3200 A— 1/^m). 

We obtained the blue channel data with the 400-line mm~^ grism, which has 
a central wave-length of 3400 A and a spectral range of 4450 A. We obtained the 
red channel data with the 400- line mm~^ grating, which has a central wavelength 
of 8500 A and a spectral range of 3800 A. Combining the blue- and red-channel 
data from the two dichroic settings, sources typically had final spectra which spanned 
the entire ~ 3200 A— 1/im optical window, albeit with higher noise at the short 
and long wavelength extremes. Based on analysis of sky lines, sources filling the 
1.5" wide slitlets used for these observations have resolution A/ A A ^ 500 and 
~ 650 for the blue and red channels, respectively. Standard stars from [43] were 
observed with the same instrument conguration for the purposes of spectrophotometric 
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Galaxy Cluster RA (J2000) Dec (J2000) z UT Date 
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calibration. Observations were generally obtained with two dithered exposures per 
dichroic conguration with typical integration times of 900 s to 1800 s, depending on 
the cluster redshift and observing conditions. This allowed both improved cosmic ray 
rejection and, by pair- wise subtraction of the images, removal of the fringing which 
strongly affects the long wavelength (A > 7200 A) LRIS data. This required minimum 
slitlet lengths of approximately 10". Since LRIS has an atmospheric dispersion 
corrector, mask position angles were optimized based on the cluster orientation and 
no special attention was necessary to align the masks with the parallactic angle. We 
show in Fig. [5] an example of our observed spectra. See the companion paper .9 for 
more details of the observations, reductions, as well as a catalog of all ~ 500 redshifts 
measured by this project. 

4.2. SDSS, SPICES and VVDS spectra 

We have extended our Keck observations with the following datasets. We take 
advantage of the improvements in calibration provided by the SDSS spectroscopic 
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Figure 5. Example of spectra from our observing program. Observed Keck 
spectra of red galaxies in the cluster MACS J2129. 4-0741 at ^ = 0.570. 



pipeline over our previous study ^221 , which was based on the early data release (EDR) 
sample, and we take advantage of the increased number of luminous red galaxies 
(LRGs) available in the public Data Release 6 (DR6) due to the larger survey area. In 
addition, we supplement our analysis using publicly available, high quality data from 
several surveys that targeted LRGs at higher redshift than the SDSS. The SPICES 
sample 04] targeted infrared-selected galaxies {K < 20) in about 1/3 square degree 
of the sky. Most of the spectra were taken at the Keck telescope. The VVDS [45 is 
the VIMOS-VLT survey carried out by the VLT/ESO telescope that spectroscopically 
targeted galaxies in the range 17.5 < Iab < 24. Neither of these surveys was aimed 
purely at LRGs, and so their efhciency at providing passively-evolving galaxies is 
significantly smaller than our Keck-survey project presented here. 
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5. Measurement of the Age-Redshift Relation 

We now proceed to measure dt/dz as a function of redshift. We fit single stellar- 
population spectra from the BC03 library for the full range of age and metallicity. We 
choose to use the BC03 library because it is the most widely used and eases the task of 
other researchers to reproduce our results. Fig.[5]shows an example of a model fit to an 
observed spectrum of a galaxy on the red envelope in the cluster MACS J0429.6— 0253. 
However, significant improvements in modeling stellar populations have taken placed 
between the release of BC03 and now, with new stellar libraries, new models and new 
differential methods that have improved the fit of the models to the data. This is an 
active field where researchers are converging to a "standard" model; until such a model 
is uncovered, the standard remains BC03. However, as we show below, the choice of 
stellar population model does not change the differential age measurements for the 
red envelope. We generate single stellar population (SSP) models with BC03 for ages 
of 0.1 — 15 Gyr and metallicities of 0.1 — 2.5 Zq. To obtain the best age, we perform 
a fit and marginalize over the metallicity to obtain error bars. Star formation is 
assumed to occur in a single burst of duration much shorter (< 5%) than the current 
age of the galaxy. We discard fits that have a reduced-x^ value > 1.2, as experience 
has shown that these galaxies have data issues (e.g., poor S/N ratio due to a bright 
moon and/or poor weather, poor sky subtraction, contamination from nearby — and, 
in four cases, lensed — sources, contamination from an active nucleus), or are better 




Figure 6. Example of a the model fit (red line) to observed spectra (black line) 
on the red envelope for a galaxy in the cluster MACS J0429. 6-0253 at z = 0.400. 
The best fit model has solar metallicity and an age of 9.6 Gyr. 
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Figure 7. Ago-redshift relation for (a) SDSS-LRG galaxies and Keck clusters; 
(b) SPICES (diamonds) and VVDS (triangles) galaxies; (c) Keck clusters; and (d) 
all samples, including the ones in Ref. 1231 . The determination of H{z) depends 
only on the determination of the upper (red) envelope; younger galaxies do not 
play any role. Note that there is a clear age-redshift relation for all samples. The 
solid line is the theoretical age redshift relation for the ACDM model. 
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Figure 8. Ages obtained witli tlie BC03 models (diamonds) and the CB08 models 
(triangles). Tlie plot shows that the red-envelope in the age-redshift relation is 
not compromised by the choice of stellar population model. The solid line is the 
theoretical age redshift relation for the ACDM model. 



fit by an extended star-formation history and therefore do not belong to the category 
of passively evolving stellar populations. We also allow for dust in the fit following 
Ref. [30] and discard those galaxies that contain dust in their fits; i.e., those where the 
dust parameter used in Ref. I30j is nonzero. To compute the differential ages, we use 
the solar metallicity, 5 Gyr-old model as our reference and compute dt with respect 
to it. 

The results of the age determination are presented in Fig. [71 where the four panels 
show the age-redshift relation for SDSS-LRG galaxies and Keck clusters (upper panel), 
SPICES (diamonds) and VVDS (triangles) galaxies (middle- upper panel), just Keck 
clusters (middle- lower panel), and all samples including the ones in Ref. [53] (bottom 
panel). Typical errors on the age are shown for in the middle- upper panel to avoid 
cluttering the plots. We also show in Fig. [H] a comparison of the age-redshift relation 
for the Keck cluster sample using BC03 (triangles) and CB08 (diamonds) models. 
Besides two discrepant galaxies, which in any case are consistent within the errors 
with the red envelope edge, the defined edge is nearly identical. This illustrates that 
the edge does not overly depend on the synthetic stellar population model. 

Note that there is a clear age-redshift relation for all samples, either taken 
individually or together. For comparison, the solid line shows the theoretical age 
redshift relation for the standard ACDM model. Also note the excellent performance 
of the clusters observed with Keck at selecting the oldest stellar populations, effectively 
from just two nights of data in good conditions. The ages derived from the Keck 
clusters nicely fill in the gap in the redshift range 0.3 < z < 1.0. Finally, we note that 
the different data samples are consistent with each other. 

We are interested in computing the "edge" of this distribution in the age-redshift 
plane. As demonstrated in Ref. [46], this is most efficiently done by computing where 
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Figure 9. The H(z) measurement from our sample of passively evolving galaxies. 
The points at z > 1 are taken from Ref. I23| . while the point at 2 = is from the 
Hubble Key Project [Sj. The dashed line is the ACDM model while the dotted 
line is an Einstein-de-Sitter model. 



the distribution, as fitted by a one-sided Gaussian, drops by 50%. The new data are 
aU at z < 1. While the recovered age-redshift relation is in good agreement with the 
SDSS one at z < 0.3, the few additional data points do not change the determination 
used in Ref. [23j. Thus, we concentrate here on the redshift range 0.35 < z < 1.0. We 
sub-divide it in two intervals: 0.35 < z < 0.6 and 0.6 < z < 1.0. The second interval 
is wider because it is populated by fewer galaxies. Each of these intervals is divided in 
two equal bins where the edge of the age-redshift relation is found in a way similar to 
that described in Ref. [13], and we also determine H{z) following Ref. [JS]. We have 
checked that the H{z) determination is robust, within the error bars, to the choice of 
binning and the choice of galaxies that populate the edge. 

6. Constraints on Cosmological Parameters 

We now proceed to place constraints on cosmological parameters using the above 
determination of H{z) in combination with the latest results from the WMAP 5-year 
release [JS] • For some applications we also consider the latest determination of Hq from 
Ref. [17]. The H{z) data help constrain cosmological parameters beyond the standard 
ACDM model, in particular for parameters that affect the Universe expansion history 
or that have degeneracies with expansion-history parameters. We find improved 
constraints for models with arbitrary curvature and dark-energy equation-of-state 
parameter w 7^ — 1, and for models with non-standard neutrino properties. 

We have obtained Monte Carlo Markov chains for the WMAP 5-year data [481149]. 
from the LAMBDA web site (www.lambda.gsfc.gov) and we importance-sampled 
them with the H{z) data. For the non-fiat dark-energy models we run independent 
Monte Carlo Markov chains. Figs. fTOHTHl illustrate how the H{z) data help to constrain 
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Table 2. H{z) from passively evolving galaxies. 





H{z) 


z 


[km sec^i Mpc^^] 





73 ±8 


0.1 


69 ± 12 


0.17 


83 ±8 


0.27 


77± 14 


0.4 


95 ±17 


0.48 


97 ±60 


0.88 


90 ±40 


0.9 


117 ±23 


1.3 


168 ± 17 


1.43 


177 ± 18 


1.53 


140 ± 14 


1.75 


202 ± 40 



cosmological parameters beyond the standard model. 

We start by considering constraints on dark energy where the dark-energy 
equation of state is parameterized as w(z) — w + Wa{^ — a) and where we allow for 
arbitrary curvature; we call this model the open-w{z)CDM model. Fig. [TOl shows the 
improvement in the w^Wa plane over CMB+Hq constraints where the Hq constraint 
comes from the Hubble Key Project [5]. Fig. [TT] shows constraints in the w-ft\ plane, 
and constraints on curvature for the above parameterization of the equation of state 
are shown in Fig. [121 Adding the Hq determination of Ref. [47] visibly improves the 
constraints only in Fig. Illl shown as the darker filled contours. 

If we impose an equation-of-state parameter constant in redshift, then the 
constraints on the curvature improve noticeably, as seen in Fig. 1131 The wide 
transparent contours are from WMAP5 data alone (with a top hat prior on Hq < 
100 km Mpc~^); the smaller transparent contours are for WMAP5 data and the 
Ref. [47] Hq determination; filled contours also add H{z) data. 

We next turn our attention to the constraints obtained on neutrino properties. 
Figs . [HI and [T5I show the constraints that can be obtained on the number of relativistic 
species, A^rei- There is a noticeable improvement over previous constraints. We find 
iVioi = 4 ± 0.5 at the la level. Finally, we present constraints on the total mass of 
relativistic species. As shown in Fig.[T6l the best constraints we obtain using the H{z) 
data is < 0.2 eV, which is an improvement of ^ 20% over the constraints obtained 
using the new Hq measurements by Ref. |47j and a factor of two improvement over 
our previous study [50] . 

Finally, we consider dark energy that is described by a slowly rolling scalar field 
and we attempt a parametric reconstruction of its potential V{z). For this purpose we 
assume the Universe to be spatially flat and we follow the steps outlined in Ref. [23] . 
We approximate a generic potential as an expansion in Chebyshev polynomials and we 
truncate the expansion at second order (i.e., we have three dark-energy parameters). 
The quantity [iJ(z)]^ is a suitably weighted integral in redshift of V{z) (see Ref. [23]; 
Eq. (29)). We then use the H{z) data presented here, the Hq determination of 
Ref. [13, and a weak Qrn prior (fi^ — 0.27 ± 0.070. The resulting reconstructed 
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Figure 11. Constraints in the w-fl/^ plane for an open-w(z)CDM model. Solid 
contours are from the WMAP5 and Hq datasets, while the solid contours use 
H{z) information. 
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Figure 12. Constraints in the ui-£7j. plane for an open-w(z)CDM model. Solid 
contours are from the WMAP5 and Hq (HST key project; [5]) datasets, while the 
solid contours use H{z) information. 




Figure 13. Constraints on w, assuming it constant, for an opon-wCDM model. 
The solid empty contours are obtained using WMAP5 data; the solid blue contours 
add the Hq constraint from Ref. |47j : the filled contours also include H(z) 
information. 
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Figure 15. Constraint on the number of relativistic species using WMAP5 
(dotted); WMAP5 with the Hq measurement of Ref. [47] (dashed); and 
WMAP5+Ho+H{z) data (solid hne). 
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Figure 16. Constraints on tlie total mass of relativistic species, in eV, using 
the WMAP5 data only (dotted line); including the new Hq measurement from 
Ref. [47) (dashed line); and adding the H{z) data which improves the Hq 
constraints by ~ 20%. The total mass of relativistic species is constraint to 
be < 0.2 eV at la. 



V{z) is shown in Fig. [iTl The reconstruction is fuUy consistent with a cosmological 
constant. 

7. Conclusions 

We report on a new measurement of the expansion history obtained from the ages 
of passively-evolving galaxies in galaxy clusters at z < 1.0. Such observations 
provide a promising new cosmological constraint, particularly for understanding the 
evolution of the dark-energy density, with a relatively modest allocation of telescope 
time, effectively just two good nights at Keck. The current measurements already 
provide valuable constraints, and the success of this first campaign should motivate 
further measurements along these lines as well as a more intensive investigation of 
the theoretical underpinnings of the calculations and the associated uncertainties. 
Further, there has been significant advancement in the last few years in modeling 
stellar populations of LRG galaxies [52l [Ml [H] and the differential technique has 
been recently applied very successfully to determine the metallicity of LRGs [S3] . It 
may be that with additional effort on both the theoretical and observational side 
that the differential-age technique may ultimately provide an important new dark- 
energy avenue which complements supernova searches, weak lensing, baryon acoustic 
oscillations, and cluster abundances. 
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Figure 17. Reconstructed dark-energy potential V{z) using an expansion in 
Chebyshev polynomials truncated at second order (i.e., three coefficients). To 
produce these constraints, we have used the H{z) data presented here, the Hg 
determination of Ref. 1471 . and a weak Qrn prior (Qm = 0.27±0.07). Gray regions 
show the 1(T and 2a confidence regions; the solid line is the best fit V{z). The 
diamonds show the redshifts at which H{z) has been measured. 
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